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Abstract. We present the first algorithm for constructing abstract 
Voronoi diagrams from bisectors that are unbounded or closed Jordan 
curves. It runs in expected O (s 2 n log(max{s, n}) Ym=2 m il^) many steps 
and 0(Y2"_ 3 mi) space, where rrn denotes the average number of faces 
per Voronoi region in any diagram of a subset of i sites. 
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1 Introduction 

Abstract Voronoi diagrams are a unifying concept that provides structural results 
and efficient algorithms for all concrete Voronoi diagrams whose bisector systems 
satisfy some simple combinatorial properties. For example, if all bisectors are 
unbounded Jordan curves, if all Voronoi regions of three sites are pathwise 
connected and if their closures cover the plane, then it is known that the Voronoi 
diagram of n sites is a planar graph of linear complexity which can be constructed 
in expected O(nlogn) many steps [4, 17, 15]. Quite a few new types of Voronoi 
diagrams were found to be under this roof, freeing authors from developing 
tailor-made algorithms [15]. 

The concept of abstract Voronoi diagrams has been extended to higher orders 
in [5,8]. Linear algorithms for special bisector systems were presented in [16, 7]. 
Also, the above conditions on the bisecting curves were relaxed in [6], by allowing 
Voronoi regions to be disconnected. But so far, only unbounded bisecting curves 
were studied, leaving many interesting Voronoi diagrams uncovered. 

In this paper we consider, for the first time, abstract Voronoi diagrams built 
from arbitrary Jordan curves, unbounded or closed. We require the following 
properties. Bisecting curves intersect vertical lines only constantly often. Any two 
bisecting curves intersect finitely often, and at most < s times if both bisectors 
are related to the same site. The closures of the Voronoi regions of any three 
sites cover the plane. Under these assumptions we can construct the abstract 
Voronoi diagram in expected 
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many steps and in 0(^™ =3 rm) space, where m, denotes the average number of 
faces per Voronoi region in any diagram of a subset of i sites. Let us look at some 
applications. 

If the bisecting curves are pseudo-circles, Voronoi regions are of linear complex- 
ity, as we will see in Section 5. Thus, < i, and (1) yields the first 0(n 2 logn) 
randomized algorithm for this class of Voronoi diagrams, which is nearly optimal 
since their size can be quadratic. 

An important special case are points with multiplicative weights; here bisectors 
are lines or Euclidean circles. The deterministic algorithm by Aurcnhammer and 
Edelsbrunner [3] runs in time 0(n 2 ) which is optimal in the worst-case where 
the diagram is of quadratic size. In this case, our algorithm takes a logn factor 
longer. But if, e.g., points are fixed, and weights drawn independently from some 
distribution (or just permuted at random) then the multiplicatively weighted 
Voronoi diagram has expected size only 0(n log 2 n), as was recently shown by 
Har-Pcled and Raichcl [14], and it can be constructed in this time. Our algorithm 
is sensitive to the smaller output size: Since m t e 0(log 2 i) holds in this case, it 
runs in expected time 0(nlog 4 n). 

For the farthest Voronoi diagram of n simple polygons of constant complexity, 
formula (1) yields a randomized (3(nlog 2 n) algorithm. Indeed, s is constant, and 
the Voronoi diagram is of linear size, as was shown by Cheong ct al. [10]; thus, 
the average number of faces per Voronoi region is constant. In comparison, the 
algorithm in [10] runs on polygons of total complexity n in deterministic time 
0{n log 3 n). 

Our algorithm for abstract Voronoi diagrams is fairly general. It uses ran- 
domized incremental construction [12,11] based on a new history graph and 
search procedure that is quite different from [17,15] and [6]. Further, it uses 
an application of Chazelles linear triangulation algorithm [9] to compute the 
trapezoidal decomposition of a given face. At a later stage this may be replaced 
by a simpler randomized method, e. g. by Amato et. al. [2], which though would 
require a new randomized analysis. 

We will discuss the new history graph and search procedure in Section 3, and 
upper bound the expected number of conflicts, and the effort invested in finding 
them in Section 4. Because of space limitations some proofs and comments are 
put in the appendix. 

2 Preliminaries 

This section contains basic definitions and well-known facts about abstract 
Voronoi diagrams. We are given a set S of n virtual sites p, q,r, . . ., and for any 
two sites p 7^ q in S a bisecting curve J(p, q) = J{q,p). These curves are supposed 
to satisfy the following five requirements. 

(Al) Under stereographic projection, each curve J(p,q) is a closed Jordan 
curve on the sphere. 
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Curve J(p, q) is unbounded iff its projection runs (in the limit) through the north 
pole of the sphere; otherwise, J(p, q) is a closed Jordan curve in the plane. Each 
such curve, unbounded or closed, bisects the plane into two open domains. As 
part of the definition of J(p, q), they are denoted by D(p, q) and D(q,p), and in 
drawings we write the site p on the D(p, q)-side of the bisector. 

Intuitively, points in D(p,q) are closer to virtual site p than to q. We define 
the Voronoi region of p e S by 

VR(p,S) := p| D(p,q), 
qes\{ P } 

and the abstract Voronoi diagram, V(S) of S, as the plane minus the union of all 
Voronoi regions. 

(A2) For any three sites pi G S, where i = 1,2, 3, the closures of the Voronoi 
regions VR(pi 1 {pi,p2,p?,}) cover the plane. 

A useful consequence of property (A2) is the following Transitivity Lemma [15]. 

Lemma 1. For any p, q,r in S we have D(p, q) f] D(q, r) C D{p, r). 

This fact would be obvious if we could read "z £ D(p, g)" as u z is closer to p than 
to </". Here it has to be derived from axioms (Al) and (A2). Lemma 1 implies 
that for every subset T C S the closures of the Voronoi regions VR(i, T) cover 
the plane [15]. 

To identify the graph structure of the point set V(S), the following property 
is convenient. 

(A3) The intersection of any two bisecting curves J(p,q),J(r,s) has a finite 
number of connected components. 

Property (A3) allows us to zoom in on any point v situated on one or more 
bisecting curves, and find a neighborhood U of v where the curves passing 
through v either stay disjoint or coincide, on either side of v. As a consequence 
one obtains a "piece-of-pie" Lemma stating that each point v of V(S) is either an 
interior point of a Voronoi edge, if v lies on the common boundary of exactly two 
Voronoi regions of different sites, or neighborhood U contains pieces of Voronoi 
regions of three different sites or more, making v a Voronoi vertex. An example 
of an abstract Voronoi diagram based on A1-A3 is shown in Figure 1. The main 
differences to standard abstract Voronoi diagrams [15] are as follows. There can 
be islands inside a Voronoi region. The boundary of an island need not contain 
a Voronoi vertex and it equals a closed bisector J(p,q). Voronoi regions can be 
disconnected. In the neighborhood of a Voronoi vertex, the same face of a Voronoi 
region can be represented several times. In order to bound the complexity of the 
Voronoi diagram we require the following. 
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Fig. 1: An abstract Voronoi diagram based on axioms Al, A2, and A3 and its 
trapezoidal decomposition. 



(A4) The intersection of any two bisecting curves J(p,q),J(p,s) related to 
the same site p has at most s connected components. 

Property A4 implies that, in a Voronoi diagram of three sites, the number of 
vertices and faces of a Voronoi region is in O(s). For the overall complexity of 
the Voronoi diagram, one obtains from (A4) the following. 

Lemma 2. A single Voronoi region is of complexity 0(sn 2 ), and this bound can 
be attained. The whole Voronoi diagram V(S) is of size 0(sn 3 ) [6]. 

The proof is by vertex counting. We observe that islands can be nested, but their 
total number is in 0(n). Namely, if an island is situated in a face of the Voronoi 
region of p, and contains points of the region of some site q ^ p, then the whole 
Voronoi region of q is confined to this island by the bisector J(q,p); see Figure 2a. 



Fig. 2: (a) The Voronoi region of q is confined to the island by the dashed bisector 
J(q,p). (b) The description of trapezoid A. 



Since we are going to subdivide Voronoi regions into trapezoids [18] we require 
one last property which is quite common in arrangement theory [19]. 

(A5) Each bisecting curve has at most a constant number of points of vertical 
tangency. All such points have pairwise different x- coordinates. 




(a) 



(b) 
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Property A5 implies that each bisecting curve has at most a constant number of 
points of vertical tangency. The trapezoidal decomposition V*(R), for a subset 
R of S, results from shooting vertical rays from all vertices and points of vertical 
tangency in both directions, until they hit another point of V(R); see Figure 1. 
Note that V*(R) is still of the same size as V(R). 

Voronoi edges and trapezoids of V* (R) are referred to by unique descriptions. 
A general edge description is shown in the leftmost drawing of Figure 3. Associated 
with vertex v is a 5-tuple containing the names of the two sites p, q separated by 
the edge, the names of the adjacent sites, r q and r p , and the coordinates of v. 
The sites are listed in clockwise order around v. A similar 5-tuple is associated 
with w. 




Fig. 3: Edges e with descriptions D R {e) = {(r q ,q,p,r p ,coov(v)), 
(u p ,p,q,u q ,coor(w))}, D R (e) = {(r q ,p,q,r p ),coor(v)} and D R (e) = {(p,q)}- 

Figure 2b depicts a general trapezoid. Its description consists of the descrip- 
tions of the upper and lower Voronoi edges e\ and ei , the sites q and I whose 
bisectors J(q,t) and J(l,t) determine the tangency points v and w, and the 
coordinates of x\ to X4. 

As a basic operation we assume that the intersections of two related bisectors 
can be computed in time O(s), see (A4). The intersection between a Voronoi 
edge of V(R) and the region of a new site VR(t, R U {t}) for t £ S \ R can also 
be computed in time O(s), because this can be decided in the diagram of 7 sites, 
the site t and the sites contained in the description of e, and e can be intersected 
in at most O(s) components. The intersection of a vertical line and a bisector 
can be computed in constant time, see (A5), and hence, the intersection of a 
trapezoid and a new region can be computed in time O(s). Further, we assume 
that for a bisector J(p, q) and a point x in the plane, we can in constant time 
decide if x is on J(p,q) or in D(p 1 q) or D(q,p). Also, for two points x and y 
on J(p, q), given a direction of J(p, q), we can decide in constant time if x or y 
comes first. 

3 Searching for Intersections 

For the algorithm we restrict our attention to the finite part of the Voronoi 
diagram and assume that a large closed curve r is given, encircling the diagram 
such that it contains all closed bisectors, it intersects each unbounded bisecting 
curve exactly twice, and no pair of bisecting curves intersect on or outside of r. 

Our algorithm is a randomized incremental construction, where the insertion 
order {n, . . . , r n } = S of the sites is given randomized, and we insert the Voronoi 
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regions successively to the current Voronoi diagram. Let R — {n, . . . ,rj}, V*(R) 
the trapezoidal decomposition of the current diagram V(R), and t e S\R the 
next site to be inserted. In order to incrementally construct the augmented 
Voronoi diagram V*(RU {t}) from V*(R) we need to determine which parts of 
V*(R) — edge or trapezoid — are intersected by the Voronoi region VR(t, R U {t}) 
of the new site t (if this region is a small island it could well be contained in the 
interior of a single trapezoid of V*(R), without intersecting any edges). 

We define edges, trapezoids, and Voronoi regions to be open, i.e., an edge 
does not contain its endpoints, and trapezoids and regions do not contain their 
boundaries. Thus, all intersections are "proper", in that they contain a subset of 
full dimension. If q were the last site inserted to the Voronoi diagram shown in 
Figure 1, its region would not intersect the previous Voronoi diagram, but one of 
its trapezoids. 

Definition 1. We say that a trapezoid A or an edge e ofV*(R) is 

(i) intersected by t, if VR(t,RU {t}) intersects A or e respectively; 
(ii) in conflict with t, if A or e with their descriptions Dr{A) and Dn(e) do not 
appear in V*(RU {t}). 

Observe that if a trapezoid or an edge is intersected by t, then it is also in 
conflict with t but not necessarily vice versa, see Figure 4a. Recall also the very 
important fact that if R' C R, then VR(i, R U {t}) C VR(t,R' U {t}), which 
follows quite directly from the definition of Voronoi regions. 




(a) (b) 

Fig. 4: (a) All trapezoids A to F are in conflict with t, because t intersects edge 
e, but only the trapezoids B to E are intersected by t. (b) Trapezoid B borders 
on two different edges of the face F. 



Finding all proper intersections of VR(i, RU {t}) with V*(R) is facilitated 
by a history graph a DAG with a single source, a vertex for each face of 

the inserted regions, and a vertex for each trapezoid ever constructed during the 
incremental process. More precisely, let R = Rj = {r\, ...,rj} be the insertion 
order, then the vertex set of H(R) equals 

{source} U {Dr. {A) \ A is a trapezoid of V*(Ri) for i £ {2, . . . , j}} 
U {F | F is a face of VR(r„ i?;) for i G {3, . . . , j}}. 
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Note that a trapezoid that exists in a sequence of Voronoi diagrams V*(Rj) is 
represented in T-L(R) by a single vertex. Faces are represented by a vertex only 
when contained in a new region, thus, in contradistinction to trapezoids, not 
all faces ever constructed are represented by a vertex in H(R). Let a vertex of 
'H(R) be called face-vertex if it refers to a face and trapezoid-vertex if it refers 
to a trapezoid. The edges of %{R) are defined incrementally, maintaining the 
invariant that the leaves of 'H(R) correspond to the trapezoids of V* (R). Also each 
trapezoid and face of the current diagram V* (R) is linked to its corresponding 
trapezoid- and face- vertex of %{R) and vice versa. 

The top of the history graph is an artificial source and its children are the 
trapezoids of the diagram of the first two sites r\ and T2 of the insertion order. 
Now let t e S \ R be the next site to be inserted. The history graph R(RU {t}) 
is obtained by updating %{R) as follows. 

If F is a face of VR(t, R U {t}), then F is made a successor of all trapezoids 
A' of V*(R) (these trapezoids are leaves of "H(-R)) intersected by F, compare 
face F in Figure 5. 

If A is a trapezoid of V*(R U {t}), which has not yet been part of V*(R) 7 
then there are two cases: 

1. Trapezoid A is contained in an r- region for an old site r € R. Then A is 
made a successor of all trapezoids A' of V*(R) intersected by A. Compare 
trapezoids B\, D\, E\, H\, H2, H3, I\, J\ in Figure 5. 

2. Trapezoid A is contained in a face F of VR(i, R U {£})• Then A is made a 
successor of F, compare trapezoids K, L, M, N, O, P in Figure 5. 



Fig. 5: Face F C VR(t, R U {t}) is inserted in the trapezoidal decomposition of 
the upper diagram and the history graph is updated. 

In addition to the history graph which contains only the descriptions of 
the trapezoids ever constructed and the "names" of the new faces, for each 
face- vertex F of R(R) we store its geometric structure. More specific, let F be 
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a face of VR(t, R\J {t}). We store its trapezoidal decomposition F* together 
with a search- structure SS(F) which allows us to perform fast point location 
for a given query point x of the plane. We link the face- vertex F of H(R) to 
its geometric structures F* and SS(F), and the trapezoids of F* and SS(F) to 
their corresponding trapczoid-vertices, which are successors of F. Further, for 
each trapezoid A of V* (R U {t}) bordering on an edge e of the boundary of F 
(this trapezoid must have been constructed while inserting F to V*(R)), we link 
the trapezoid- vertex A to the edge e in F* , and vice versa, compare trapezoids 
Ei,Hi,H2,Hs,I\, Ji in Figure 5. Note that trapezoid A may border on up to 
two different edges of the boundary of F, compare Figure 4b. If this is the case, 
we link A to both of them. 

Now let us discuss how to compute the geometric structures of F from scratch. 
Each edge on the boundary of F can be split up into a constant number of x- 
monotone arcs, by putting an additional vertex on each point of vertical tangency 
(A5). By doing this, in time 0(|-F|), we obtain a simple polygon with monotone 
curved edges whose trapezoidal decomposition can be computed in time 0(|F|), 
see [9]. Observe that we consider F without the rest of the diagram and we 
compute the trapezoidal decomposition both in the inside and outside of F. The 
trapezoids in the outside are needed to obtain a monotone subdivision F* of 
the whole plane. For this subdivision we can now built a search-structure SS(F), 
based on the structure for monotone subdivisions introduced by Edelsbrunner 
et. al. in [13], which again takes time and storage 0(|F|), and it allows us to 
perform point location for a given query point x of the plane in time 0(log \F\). 

Both in [9] and [13] the edges are straight lines, but the same ideas can be 
adapted to instances where edges are monotone curves, and the intersection 
between an edge and a vertical line can be computed in constant time. 

Now we describe how to find all trapezoids of V*(R) intersected by the 
region of the new site t. We walk through the history graph along its trapezoid- 
vertices intersected by t as follows, starting with the successors of the source. If 
a trapezoid- vertex A of V*(Ri) is intersected by VR(i, Ri U {t}), then we test 
recursively all its succeeding trapezoid-vertices for intersection with t; for the 
succeeding face-vertices we discuss later how to proceed. First, if A is linked to 
an edge e of the trapezoidal decomposition F* of a face- vertex F, we test if e is 
intersected by t and if yes, we test all trapezoids bordering on e from the inside 
of F for intersection with t; these trapezoids are found in F* . Afterwards we set 
a flag on e saying that e already has been tested for intersection with t, thus if e 
is reached from another trapezoid, we do not test it and its adjacent trapezoid for 
intersection with t again. For the running time analysis later on we note already 
now that if edge e is intersected by t, then all trapezoids bordering on e (there 
may be many!) are in conflict with t. 

For each trapezoid bordering on e from the inside of F, having been successfully 
tested for intersection with t, we recursively test also its at most 4 neighbors in F 
for intersection with t. Doing this we trace the region of t through F. Of course 
the region of t may intersect F and other faces of the same region in several 
disconnected components, but we will later see that we will find a starting point 
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Fig. 6: The region of t (faces with normal boundary) is traced through the 
trapezoidal decomposition F* (face with fat boundary), in the left figure from 
two starting points - one from A trough e\, and one from B through and in 
the right figure from a point x found using SS(F). 

for tracing each component. Again we set a flag with respect to t on each tested 
trapezoid to prevent it from having to be tested again. 

In the left drawing of Figure 6, the region of t intersects F* in 3 components, 
F\,F2, and F3. If we enter F* from A through e\, then all trapezoids of F* 
intersected by F\ and F 2 are found. To detect the trapezoids intersected by F 3 , 
we need another starting point, e.g. from B through e 2 . 

We may now have detected additional trapezoids which are intersected by 
t, they are linked to their corresponding trapezoid in the history graph and we 
recursively test also their successors for intersection with t. 

For a face-vertex F r . C VR(rj,.Rj) reached in the history graph, we choose 
an arbitrary point x of the bisector J(t,ri) and use SS(F ri ) to determine if x is 
contained in F and if yes, in which trapezoid. If no trapezoid of F ri containing x 
is detected, i. e. x lies outside of F ri , then nothing has to be done, as will become 
clear later. Otherwise a trapezoid A C F r . is detected with x £ A. Then we know 
that A is intersected by t, and like before we can trace the region of t through F r . 
and mark the tested trapezoids with flags with respect to t (unless they already 
have a flag with respect to t, then we are done), compare the right drawing of 
Figure 6. Again we may have detected additional trapezoids intersected by t and 
we recursively test also their successors in the history graph for intersection. 

Lemma 3. By walking through the history graph H(R) as described above, we 
reach all trapezoids of V* (R) which are intersected by VR(t,RU {t}). 

Proof. (See first proof in Appendix) We prove that if a trapezoid A of H(R) is 
intersected by t, then either a predecessor of A in H(R) is intersected by t, or A 
is detected in the search structure SS(F) of a face F, or found while tracing the 
region of t through a face F. In all these cases, we will detect A in the history 
graph. □ 

Let E t be the set of all trapezoids of V*(R) intersected by t. Once this set 
has been computed, we can update V*(R) like described in [6] and %{R) like 
described above. This takes time 0(s 2 \E t \), see [6]. 
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4 Analysis 



To organize the discussion of the running time of our algorithm we define the 
following variables. Like before let Rj = {n, . . . , rj} C S, the sites are inserted 
in this order, and t = r J+1 e S \ Rj is the next site to be inserted. Further 

E t := {A | A is a trapezoid of V*(R) intersected by t}, and 

j 

c := ^ # trapezoids of V*(Ri) not in in conflict with t. 

i=3 

Lemma 4. The set E t can be computed in time 0(s 2 c log(max{s, n})). 

For the following randomized analysis we use the techniques by Clarkson, 
Mchlhorn, Seidel [11]. The proofs are similar to the ones in [6], for completion 
the missing ones can be found in the Appendix. To upper bound the number 
of conflicts during the insertion process in the algorithm, we define mi as the 
average number of faces per region, over all diagrams of i sites from S. 

Lemma 5. For R = {n, ...,rj}, the expectation of c is 0(J^_ 2 ^). 

Lemma 6. The expected size of the history graph 'H(-S') is 0(Y^ =3 mi). 

Lemma 7. The expected size of E t is 0(mj). 

Theorem 1. Let {J(p,q) : p ^ q € £} be a bisecting curve system fulfilling 
axioms (Al) to (A5). Then V{S) can be computed in expected time 

n 

o(s 2 nlog(max{s,n})^— ), (2) 

i=2 1 

and expected space 0(^" =3 mj). 

Proof. Because of Lemma 4 and 5, the sets E Ti , 3 < i < n can be computed in the 
claimed expected time, formula (2). Afterwards H(S) and V*(S) can be computed 
in time 0(s 2 J27=3 l-^nl)> see the end of the previous section, and because the 
expectation of \E r . \ is 0(rrii), the expected value is in 0(s 2 ^" =2 mi) e (2). 

Lemma 6 shows the expected size of H(S), and the size of the trapezoidal 
decomposition and the search structures for the faces F is linear in the size of F, 
see [13], which gives us the expected space. □ 

5 Pseudo-circles: an Application 

Let the set of bisectors J :— { J(p, q) : p ^ q e S} be a set of pseudocircles. 
Thus, each bisector is a simply closed curve and any two bisectors J(p, q) and 
J(r, t) are either equal, or their intersection is empty, exactly one nontransversal 
point, or exactly two transversal points. Such a curve system fulfills axioms (Al), 
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(A3), and (A4), in addition we assume that also (A2) and (A5) are fulfilled. 
If | S\ = 3 then all Voronoi regions are connected, but in contradistinction to 
Voronoi diagrams with only unbounded bisectors, this does not imply that all 
regions in V(S) are connected, see the Appendix for more details. 

Nevertheless, we have the following complexity bound on the size of the 
Voronoi diagram. Agarwal et. al. [1] showed that the boundary of the union of 
n pseudodiscs, for n > 3, consist of at most 6n — 12 elementary arcs. In our 
case we are interested in the complexity of the intersection of pseudodiscs and 
complements of pseudodiscs, thus the result by Agarwal et. al. is not directly 
applicable, but the technique used in the last part of our proof is similar. 

Theorem 2. Let S be a set of n sites together with a set of bisecting curves 
fulfilling our axioms. Then each Voronoi region in V(S) is of size 0(n) and V(S) 
is of size 0(n 2 ). These bounds are tight in the worst case. 

Proof. (Sketch) Let p be a site in S. Its Voronoi region is the intersection of open 
pseudodiscs and complements of pseudodiscs. 

Claim 1. The intersection of n open pseudodiscs is simply connected 1 (or empty). 
Claim 2. If each edge has the name of the pseudocircle which it is part of, 
then the edges along the boundary of the intersection of n pseudodiscs is a 
Davenport-Schinzel-Sequence of order 2. 

Claim 3. The intersection of n — 1 open pseudodiscs and complements of pseu- 
dodiscs has O(n) edges on its boundary. □ 

Theorem 3. The Voronoi diagram based on a bisector system of pseudocircles 
fulfilling axioms (A2) and (A5) can be computed in expected time 0(n 2 log n). 

6 Conclusion 

Our algorithm covers many concrete Voronoi diagrams for which yet no efficient 
algorithm existed. The run time strongly depends on the size of m,j, the average 
number of components of a Voronoi region in a diagram of j sites. As we have 
seen, for pseudocircles this works very nicely and we obtain a nearly optimal run 
time, but there are examples were nij is higher than the total complexity of the 
final diagram, thus in these cases our algorithm would be less efficient. 

Also, if one wants to implement the algorithm one may not want to use 
Chazelles linear triangulation algorithm [9], but an easier approach, e.g. a 
randomized one [2] . This may be possible but it would require a new randomized 
analysis of the run time. 
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7 Appendix 



7.1 Proofs of Section 3 

Proof, (of Lemma 3) We prove that if a trapezoid A of H(R) is intersected by t, 
then either a predecessor of A in H(R) is intersected by t, or A is detected in 
the search structure SS(F) of a face F, or found while tracing the region of t 
through a face F. In all these cases, we will detect A in the history graph. 

So, assume that A is a trapezoid which was constructed during the insertion of 
a site Ti G R, i. e. A is in V*(Ri) but not in V* (Ri\{ri}). Now there are two cases. 

Case 1: A C VR(r, RA and r^r,. 

Let Ai,...,Ak be the predecessors of A in H(Ri). Then A C U ie ^ fc j Aj. 
Thus if A is intersected by t, then also at least one of its predecessors must be 
intersected by t. 

Case 2: A C VR(r;,i?;). 

Let F be the face of VR(rj, i£j) containing A. 

If the region of i is not completely contained in F, then no face of the region 
of t may be an island in F. Thus there is a corridor contained in the region of 
t starting in A , crossing an edge e on the boundary of F and reaching into an 
adjacent trapezoid B. The trapezoid B must have been constructed during the 
insertion of rj, thus the intersection of t and B is detected like in Case 1. Because 
the boundaries of B and F share a common edge e which is intersected by t, the 
region of t is traced through F* and A is detected. 

Now let the region of t be completely contained in F. Lets first observe 
that F is a face-vertex in H(Ri) and if A\, . . . , A^ are its predecessors, then 
F is contained in Uie{i k} ^i- Thus at least one of the predecessors of F is 
intersected by t and we reach F together with its search structure SS(F) while 
walking through H(R). Furthermore, because the region of t is contained in F 
the whole bisector J(r i; t) must be contained in F, thus for any point x on the 
bisector we will detect a trapezoid B contained in F, which is intersected by t. 
From there we start to trace the boundary of the region of t and reach A. □ 

7.2 Proofs of Section 4 

Proof, (of Lemma 4) To test if a trapezoid A or edge e is intersected by t takes 
time O(s), using our basic operations. Each trapezoid-vertex of "H(i?) has an 
outdegree of O(s), because each trapezoid can be intersected in at most O(s) 
components by a new Voronoi region. Thus for each intersected trapezoid-vertex 
at most O(s) succeeding trapezoids are unsuccessfully tested for intersection with 
t. This takes time 0(s 2 ) per intersected trapezoid vertex. 

Further, each trapezoid-vertex is linked to at most two edges e of the trape- 
zoidal decomposition F* for a face-vertex F. To test them for intersection with t 
or to detect that they already have a flag with respect to t takes time O(s) per 
intersected trapezoid-vertex. 
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If such an edge e has been successfully tested for intersection with t, we have 
to test all trapezoids bordering on e from inside of F for intersection with t. 
There may be many such trapezoids, but fortunately, if e is intersected by t, then 
all trapezoids bordering on e are in conflict with t and thus they can be upper 
bounded by c. Within the face F, each trapezoid has at most 4 neighbors, thus 
while continuing the tracing within the face F, for each intersected trapezoid at 
most four adjacent trapezoids may be intersected unsuccessfully for intersection 
with t. The flags prevent us from testing edges and trapezoids several times for 
intersection with t. Thus altogether this process takes time O(sc). 

For each face- vertex F of H(R) reached, we have to perform a point loca- 
tion searching in SS(F) which takes time 0(log(F)). A very rough (but suffi- 
cient for O-notations) estimation of \F\ is 0(sn 2 ), which gives us 0(log(F)) e 
O(log(max{s, n})). If the searching was successful, we trace the region of t and 
again for each intersected trapezoid we test at most 4 additional trapezoids 
unsuccessfully. 

Each trapezoid intersected by t is also in conflict with t, thus the overall 
running time to compute E t is in (9(s 2 clog(max{s, n})). □ 

Proof, (of Lemma 5) Let A be a trapezoid- vertex of the history graph H(ri , . . . , r j ) , 
where the sites n, . . . , rj are inserted in this order. Now let t — r J+ i e S \ R 
randomly chosen. If A is in conflict with t, then A with its description Dr{A) 
would never have been constructed, if t would have been inserted as the first site, 
i. e., A $l H(t, n,..., rj). By elementary set calculus, it follows: 

\H(r 1 ,...,r j )\H(t,r 1 ,...,r j )\ 
=\H(n,. . . , rj)\ - \H(t, ri , ...,rj)\ + \H(t, n, . . . ,rj) \ H(n, . . .,rj)\. 

The expectation of \H{r\, . . . , rj)\ — \H(t, n, . . . , rj)\ is < 0, because |{n, . . . , rj}\ < 
\{t, n,..., rj}\ implies for randomized inputs the expectation \%{r\, . . . ,rj)\ < 
\H(t, ri, . . . , rj)\. Thus the equation above has an expected value 

<E(\H(t,r 1 ,...,r j )\H(r 1 ,...,r J )\). 
s ^ ' 



If A G \H(t, n, . . . , rj)\H(ri, . . . , rj)\, then t € set(D RU { t }(A)). Let i be minimal 
with A e V*(t,ri, . . . ,n) and thus n G sct(L> flu{t } (A)). 
Then the size of X is upper bounded by 

3 

< V \{A e V*(t, n , . . .,n)\t, n e set(D RiU{t} (A))}\, 

.■ o " ' 
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where Ri = {ri, . . . , r,}. By choosing t and ti randomly from {t, n, . . . , r^} we 
can estimate the expectation of Y by 

^T^Ttv E \{AeV*(t,n,...,n)\x,y€set(D RiU{t} (A))}\ 

(a:,j/)e{t,ri,...,ri} 2 ,x^y 

and since set(D R (A)) consists of at most 16 different sites it follows 

<7tt^v-16-15- \V*(t, ri ,...,n)\ <C^ 
[i + 1)1 - v ' i 

<C\V(t,r 1 ,...,r i )\<Cm i+1 (i+l) 

for a constant C. Thus the expectation of X is 

i=2 * \i=2 * / 

This implies c e 0(EL2 f)- D 

Proof, (of Lemma 6) Let i?i = {n, . . . , r^} be the set of the first i sites inserted in 
this order. The size otH(S) equals £™ =3 H(Ri)\H(Ri-i) . If a; G H(Ri)\H(Ri-i) , 
then by the definition of the history graph, either x is a a face of VR(r i7 ors 
is a trapezoid of \ V*(Ri-i). Because the number of faces in the history 

graph is less than the number of trapezoids (each face-vertex has at least one 
unique succeeding trapezoid) , the number of all vertices of the history graph is 
< two times the number of trapezoid-vertices. If a; is a trapezoid-vertex, then 
e D R . (x) which implies: 

s(h(^)\^-i))<t E E 5 ( A ' r )' 

reRi .4etrapezoids(V'*(i?i)) 



where 5(A, r) 



1, if resets (A)) 
0, else. 

<?(16|V*(i2i)|) e 0(mi). 



Summing up over all i from i = 3 to n, we get E(\H(S)\) £ 0(^" =3 mi). □ 

Proof, (of Lemma 7) The set E t equals the set of all trapezoids of V*(Rj) \ 
V*(Rj U {t}). With the same arguments as in the two previous proofs one can 
show that the expected size of this set is 0(rrij). □ 



7.3 Comments about Section 5 

Lemma 8. If \S\ = 3, then all Voronoi regions are connected. 

Proof. Let S — {p,q,r}. Then the region of p equals the intersection of the 
two dominance regions D(p, q) n D(p,r). A dominance region is either an open 
pseudodisc or the open complement of a pseudodisc. It is easy to see that the 
intersection of two dominance regions is always connected. □ 
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For unbounded bisecting curves we have the property that if for all subsets 
5" of S of size 3 the Voronoi regions in V(S') arc connected, then all Voronoi 
regions in V(S) are connected, see [15, Lemma 14]. This is something which 
is no longer true, if bisectors may be closed curves, as has been shown already 
for multiplicatively weighted Voronoi Diagrams [3]. Here one region may consist 
of 0{n) faces and the total complexity of V(S) is 0(n 2 ). Interestingly there is 
no connection at all between the property that for all subsets S' of S of size 
k, 1 < k < n — 1, all regions in V(S') are connected and the property that all 
regions in V(S) are connected. 

Lemma 9. Let S be a set of n sites. We can choose a set of bisecting curves 
such that for all proper subsets S' C S all Voronoi regions in V(S') are connected 
but there is a region in V(S) which is disconnected. 

Proof. An example for such a system of bisecting curves is depicted in Figure 7. 
Let S = {pi, . . . ,p n }- The uppermost circle is the bisector J{p\,P2), J{P2,Pz)i 
. . . , J(p2,Pn), where D{pi,p 2 ), D(p 3 ,p2) ■ ■ ■ , D(p n ,p2) equal the outer face of the 
circle. The circle to the right of the uppermost one is the bisector J{p\ 1 pz), 
J(P3,Pi), ■ ■ -,J(P3,Pn), where D(pi,p 3 ), D(p 4 ,p 3 ) . . .,D(p n ,p 3 ) equal the outer 
face of it. Finally the circle to the left of the uppermost one is the bisector 
J(pi,p n ), where D(pi,p n ) is the outer face. 

Now the region of p\ equals the plane minus the union of the circles. The 
uppermost disc is the region of pi , the disc to the right minus its left one is the 
region of p 3 and so on. Thus in V(S), where all circles are present, the region of 
Pi is disconnected into two faces. But as soon as one circle is missing, i. e. when 
we have a subset S' not containing one of the sites p 2 ,...,p n , then the region of 
Pi grows together through this circle and becomes connected in V(S'). All other 
regions are connected in all V(S'), where S' is a subset of 5\ □ 




Fig. 7: An example where the region of pi is disconnected in V(S), but for proper 
subsets S'CS all regions in V(S') are connected. 



1G 



Proof, (of Theorem 2) 

Let p be a site in S. Its Voronoi region is the intersection of open pseudodiscs 
and complements of pseudodiscs. 

Claim 1. The intersection of n open pseudodiscs is simply connected 2 (or empty). 

Proof. It is clear that the intersection has no holes, otherwise one of the pseu- 
dodiscs would have a hole. Thus, it remains to show that it is connected, which 
follows from [15, Lemma 14], because if all Voronoi regions in any diagram of 
three sites are connected, then related unbounded bisectors may intersect in at 
most two points, and these intersections are transversal. Thus, the bisectors can 
be connected to closed ones at infinity giving us a system of pseudocircles. Now 
the claim follows from the lemma. □ 

Because of Claim 1, the boundary of the intersection of n pseudodiscs is a 
closed curve consisting of pseudocircle segments. We call these segments edges. 

Claim 2. If each edge has the name of the pseudocircle which it is part of, 
then the edges along the boundary of the intersection of n pseudodiscs is a 
Davenport-Schinzel-Sequence of order 2. 

Proof. Let an edge have the name q if it belongs to the pseudocircle J(p, q). Now 
suppose there are edges q and r alternating three times in the sequence, i. e. 
q...r ...q...r. Then the two bisectors J(p, q) and J(p, r) must intersect in more 
than two points, a contradiction. □ 

This claim shows that there can be at most 2n — 1 edges on the boundary of 
the intersection of n pseudodiscs. 

Now recall that VR(p, S) is the intersection of both open pseudodiscs and 
complements of pseudodiscs. Let F be the intersection of the pseudodiscs D(p, qj). 
By Claim 2, F is connected and has 0(n) edges on its boundary. It is clear that 
VR(p, S) is contained in F. Place a vertex on each edge on the boundary of F 
and for each pseudodisc D(q,p) (whose complement is a dominance region of p) 
intersecting F, place a vertex q in D{q 1 p) n F; observe that each pseudodisc can 
intersect F in at most one connected component. 

Connect two vertices q and r belonging to D(q,p) and D(r,p) by an edge iff 
dD(q,p) and dD(r,p) intersect in a point in F not contained in any other D(s,p), 
i. e. J(q,p) and J(r,p) intersect in a vertex on the boundary of the region of p 
Further, connected a vertex q belonging to D(q,p) and a vertex v belonging to 
an edge e on the boundary of F iff dD(q,p) intersects e in a point not contained 
in any other D(r,p). Finally, connect two vertices belonging to edges e and e' on 
the boundary of F, if e and e' are incident. 

Now consider the resulting graph G, see Figure 8. Because two pseudocircles 
can intersect in at most two points, G is a planar graph without parallel edges. 

2 We call a set simply connected if it is both connected and has no holes. 
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Fig. 8: Proof of the last part of Theorem 2. Fat curves depict the boundary of 
F, normal curves the boundaries of the pseudodiscs D{q 1 p) and dashed curves 
edges of the graph G. 

Further we can draw the edges e between two vertices q G D(q,p) and r G D(r,p) 
such that e C D(q,p) U D(r,p), which is possible, because by definition e exists 
only if D(q,p) P\D(r,p) ^ 0. Then it is clear that each connected component of 
VR(p, S) corresponds to a unique face of G (but G may have more faces than 
VR(p,S)). Because G has 0(n) vertices it also has 0(n) faces and edges, which 
proves the upper bound of the Theorem. 

An example, where the bounds are tight is the Voronoi diagram of multiplica- 
tively weighted points, see [3]. □ 
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